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A new systematic calculation of magnetization and specific 
heat contributions of vortex liquids and solids (not very close 
to the melting line) is presented. We develop an optimized 
perturbation theory for the Ginzburg - Landau description of 
thermal fluctuations effects in the vortex liquids. The expan- 
sion is convergent in contrast to the conventional high tem- 
perature expansion which is asymptotic. In the solid phase 
we calculate first two orders which are already quite accu- 
rate. The results are in good agreement with existing Monte 
Carlo simulations and experiments. Limitations of various 
nonperturbative and phenomenological approaches are noted. 
In particular we show that there is no exact intersection point 
of the magnetization curves both in 2D and 3D. 



It was clearly seen in both magnetization Q and spe- 
cific heat experiments |Q that thermal fluctuations in 
high Tc superconductors are strong enough to melt the 
vortex lattice into liquid over large portions of the phase 
diagram. The transition line between the Abrikosov vor- 
tex lattice and the liquid is located far below the mean 
field phase transition line. Between the mean field tran- 
sition line and the melting point physical quantities like 
the magnetization, conductivity and specific heat depend 
strongly on fluctuations. Several experimental observa- 
tions call for a reflned precise theory. For example, a 
striking feature of magnetization curves intersecting at 
the same point {T*,H*) was observed in a wide rage of 
magnetic fields in both the layered j|] materials and the 
more isotropic ones [Q. To develop a quantitative the- 
ory of these fluctuations, even in the case of the lowest 
Landau level (LLL) corresponding to regions of the phase 
diagram "close" to Hc2 is a very nontrivial task and 
several approaches were developed. 

Thouless and Ruggeri Q proposed a perturbative ex- 
pansion around a homogeneous (liquid) state in which 
all the "bubble" diagrams are resummed. Unfortunately 
they proved that the series are asymptotic and although 
first few terms provide accurate results at very high tem- 
peratures, the series become inapplicable for LLL di- 
mensionless temperature qt ^ [T — Tmf{H))/{TH)^/'^ 
smaller than 2 in 2D quite far above the melting line 
(believed to be located around — ~12). Generally at- 
tempts to extend the theory to lower temperatures by the 
Borel transform or Pade extrapolation were not success- 



ful [0. Several nonperturbative methods have been also 
attempted including renormalization group and the 
expansion ||^. Tesanovic and coworkers developed 
a theory based on separation of the two energy scales 
|]l0| : the condensation energy (98%) and the motion of 
the vortices (2%). The theory explains the intersection 
of the magnetization curves. 

In the first part of paper we apply optimized perturba- 
tion theory (OPT) first developed in field theory ^1 12 
to both the 2D and 3D LLL model. It allows to obtain 
a convergent (rather than asymptotic) series for magne- 
tization and specific heat of vortex liquids together with 
precision estimate. The radius of convergence is ot — ~3 
in 2D and = — 5 in 3D. On the basis of this one can 
make several definitive qualitative conclusions. 

Our starting point is the Ginzburg-Landau free energy: 
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where A = [By, 0) describes a nonfluctuating constant 



magnetic field in Landau gauge and D = V - 
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Lc is the width (for simplicity we write expressions 
for the 2D case, essential 3D complications are discussed 
separately). For simplicity we assume a(T) = aTc(l — t), 
t = T/Tc- On LLL, the model after rescaling reduces to 
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where the LLL reduced temperature ar = v JxJ 
is the only parameter in the theory Q. Here b 

We will use a version of OPT, the optimized gaussian 
series jl^ . It is based on the " principle of minimal sensi- 
tivity" idea [ pT| , first introduced in quantum mechanics. 
Generally a perturbation theory starts from dividing the 
Hamiltonian into a solvable "large" part K and a pertur- 
bation V. Since we can solve any quadratic Hamiltonian 
we have a freedom to choose "the best" such quadratic 
part. Quite generally such an optimization converts an 
asymptotic series into a convergent one (see a compre- 
hensive discussion, references and a proof in fl^ ). 

Due to the translational symmetry of the vortex liq- 
uid there is just one variational parameter, e, in the free 
energy divided as follows: 
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where an = ax ~ £■ One reads Feynman rules from 
eq.(||): K determines the propagator (just a constant), 
the first term in y is a "mass insertion" vertex with a 
value of -!-arr while the four line vertex is To cal- 
culate the effective free energy density f^ff = — 47rlnZ, 
one draws all the connected vacuum diagrams. We cal- 
culated directly diagrams up to the three loop order. 
However to take advantage of the existing long series of 
the non optimized gaussian expansion, we found a re- 
lation of the OPE to these series. Originally Thouless 
and Ruggeri calculated these series /e// to sixth order, 
but it was subsequently extended to 12*'* (9*'* in 3D) by 
Brezin et al and to 13*'' by Hu et al It is usually 

presented using variable x introduced by Thouless and 

Ruggeri [|| x = — ^ (^ar + y^a'^ + 16^ as follows: 

f^^j = 21og^ -I- 2J2^=iCnX"- We can obtain all the 
OPT diagrams which do not appear in the gaussian the- 
ory by insertions of bubbles and mass insertions from 
the diagrams contributing to the nonoptimized theory. 
Bubbles or "cacti" diagrams are effectively inserted by a 
technique known in field theory [[l4[: 



feff = 2 log ^ -f 2 ^ c„ 
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Summing up all the insertions of the mass vertex is 
achieved by £2 = e + aan- Here a was introduced to 
keep track of order of the perturbation, so that expand- 
ing feff to order a"^^, and then taking a = 1 we obtain 
/„(e) (calculating /„ that way, we checked that indeed 
the first three orders agree with the direct calculation). 
The n*^ OPT approximant /„ is obtained by minimiza- 
tion of /n(e) with respect to e: 
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The above equation is equal to 1/ £^"+3 times a polyno- 
mial g„ (z) of order n in z = e • an- That eq.(P) is of 
this type can be seen by noting that the function / de- 
pends on the combination a/ (e -\- aan)^ only. We were 
unable to prove this, but have checked it to the 40*'' or- 
der. This property greatly simplifies the task: one has 
to find roots of polynomials rather than solving transcen- 
dental equations. There are n (real or complex) solutions 
for gn (z) = 0. However (as in the case of anharmonic 
oscillator [Q) the best results gives a real root with 
the smallest absolute value. We then obtain e{aT) = 

\ {aT + v^fly — 4z„^ solving 2„ = e • an = £aT — 

On Fig. 1 we present OPT for different orders includ- 
ing n — (gaussian) together with several orders of the 



nonoptimized high temperature expansion. One observes 
that the OPT series converge above ax = —2.5 and di- 
verge below aT = —3.5. The proof of convergence is anal- 
ogous to that for the anharmonic oscillator, see ref. [^ . 
On the other hand, the nonoptimized series never con- 
verge despite the fact that above ct = 2 first few approx- 
imants provide a precise estimate consistent with OPT. 
Above ot = 3 the liquid becomes essentially a normal 
metal and fluctuations effects are negligible (see Fig. 2, 
3). Therefore the information the OPT provides is essen- 
tial to compare with experiments on magnetization and 
specific heat. If precision is defined as (/12 — /lo) //lo, 
we obtain 4.87%, 1.27%, 0.387%, 0.222%, 0.032% at ar = 
—2,— 1.5,— 1,-0. 5,0 respectively. For comparison with 
other theories and experiments on Fig. 2 and 3 we use 
the 10*'' approximant. 

The calculation is basically the same in 3D, the only 
complication being extra integrations over momenta par- 
allel to the magnetic field. However since the propaga- 
tor factorizes, these integrations can be reduced to cor- 
responding integrations in quantum mechanics of the an- 
harmonic oscillator ^,|ll|. The series converge above 
aT = —4.5 and diverge below ot = —5.5. The nonop- 
timized series are useful only above ax — — 1 . The agree- 
ment is within the expcted precision when we compare 
our results in 3D with ref. ||l5| . 

Now we turn to the vortex solids. Here the minimiza- 
tion is significantly more difficult due to reduced sym- 
metry. Unlike in the liquid the field %p acquires a non- 
homogeneous expectation value and can be expressed as 
'4>{x) = v{x) -\- x(a;), where x describes fluctuations. As- 
suming hexagonal symmetry, it should be proportional 
to the mean field solution v{x) = vipk=o{x) with a vari- 
ational parameter v taken real thanks global U{1) gauge 
symmetry where tpk (x) is the quasi - momentum basis on 
LLL m. Expanding x 
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exp[-i9k/2]ipk{x) (Ok + lAk) . (6) 



where real fields Ak = A*_^ {Ok = 0*_^) describing 
acoustic (optical) phonons of the flux lattice. The phase 
exp[— i0fe/2] defined, as in the low temperature perturba- 
tion theory developed recently , via 7^ = \^k \ exp[i0fc] , 
7fe = {ip{){x)ipo{x)(pl{x)(p*^{x)) ^, is crucial for simplifica- 
tion of the problem. The most general quadratic form 
is 

I OkGohik)0^k + AkG^\{k)A_k + 

iSTT Jk 

OkGo\{k)A^k + AkGo\ik)0^k, (7) 

with matrix of functions G{k) to be determined together 
with the constant v by the variational principle. The 
corresponding gaussian free energy f^f / is 

qtv + —V — 2 
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+ {v^ [(2/3fe + |7fe|) Goo (k) + (2/3fc - |7fe|) Gaa (k)]), 
+ iPk-i [Goo (fc) + Gaa (k)] [Goo (0 + Gaa (01)^,; 

+ ^ {(l7fc| (Goo (k) ~ Gaa {k)))l + 4 (|7,| Goa {k))l] 

where (•■•)/. denotes average over Brillouin zone (3k = 
{iPo{x)ipoix)ipl{x)(pk{x)) ^ ,/3a = Po- The gap equations 
obtained by the minimization of the free energy look 
quite intractable, however they can be simplified. The 
crucial observation is that GoA(k) = is a solution and 
general solution can be shown to differ from this simple 
one just by a global gauge transformation. One can set 
'E{k) + A\jk[ 

E{k)-A\-fk[ 
A is a constants (details will appear elsewhere). The 
function E{k) and the constant A satisfy: 
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Observing that (3^ has a very effective expansion in 
X = exp[-a^/2] = exp[-27r/V3] = 0.0265, ^fc = 
Ejr=oX"/5n(fc), /3„(fc) = X;|xp=„ai exp[ik»X] and us- 
ing the hexagonal symmetry of the spectrum, E(}i) can 
also be expanded in "modes" E{U) — J2 ^nPnik)- The 
integer n determines the distance of a points on the 
hexagonal lattice X from the origin. One estimates that 
En — x"flT, therefore the coefficients decrease exponen- 
tially with n. For some integers, for example, n = 2,5, 6, 
(3n = 0- We minimized numerically the gaussian energy 
by varying v,A and first few modes of E(k). In prac- 
tice two modes are quite enough. The results show that 
around < —5, the gaussian liquid energy is larger than 
the gaussian solid energy. So naturally when ax < —5, 
one should use the gaussian solid to set up a perturbation 
theory instead of the liquid one. The gaussian energy in 
either liquid (see line TO on Fig.l) or solid is a rigor- 
ous upper bound on the free energy. We calculated the 
leading correction (without its minimization) in order to 
determine the precision of the gaussian result (see Fig. 3 
for the specific heat results). We obtain 0.2%, 0.4% and 
2% at ar = -30, -20, -12 respectively. 

In the rest of the paper we compare our results with 
other theories, simulations and experiments. An analytic 
theory used successfully to fit the magnetization and the 
specific heat data jlTj was developed in Their free 
energy density is: 
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The corresponding magnetization and specific heat are 
shown as a dashed lines on Fig. 2 and 3 respectively. 
At large positive ar, feff = 21ogaT + ^ - ^ + |^ 
and differs very little from the exact series 2 log ay + 
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It's low temperature asymptotics is how- 



ever less precise:— 2^ — 21og^^ which has an oppo- 
site sign of the log term compared to the exact series 
|l|-^-h21og^-^. This is seen on Fig.3 quite 

clearly. Instead of rising monotonously from G/AC — 1 
till melting as is predicted by OPT, their curve (dashed) 
first drops below 1 and only later develops a maximum 
above 1. In the liquid region it underestimates the specific 
heat. We conclude therefore that although the theory of 
Tesanovic et al is very good at high temperatures they 
become of the order 5 — 10% at ut = —3. An advantage 
of this theory is that it interpolates smoothly to the solid 
and never deviates more than 10%. 

Experiments on great variety of layered high Tc 
cuprates {Bi or Tl [|j based) show that in 2D, magnetiza- 
tion curves for different applied fields intersect at a single 
point (M*,T*). The range of magnetic fields is surpris- 
ingly large (from several hundred Oe to several Tesla). 
This property fixes the scaled LLL magnetization defined 

as m(aT) = -^SsuM = in^ f^M. Demanding that 

the first two terms in l/a^ expansion of m{aT) are con- 
sistent with the exact result, one obtains 



"1 (ot) = ^ ( ot - a/16 
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When it is plotted on Fig. 2 (the dotted line), we find 
that at lower temperatures the magnetization is overes- 
timated. The OPE results are consistent with the exper- 
imental data (points) within the precision range till 
the radius of convergence aj- = —3. It is important to 
note that deviations of both the phenomenological for- 
mula eq. ([l0|) and the Tesanovic's are clearly beyond our 
error bars. Therefore we conclude that the coincidence of 
the intersection of all the lines at the same point {T* , M*) 
cannot be exact. Like in 3D the intersection is approx- 
imate, although the approximation is quite good espe- 
cially at high magnetic fields. 

Specific heat OPE result in 2D is compared on Fig. 3 
with Monte Carlo simulation of the same model by Kato 
and Nagaosa |]l^ (black circles) (arid the phenomeno- 
logical formula following from cq. (p^ , dotted line). The 
agreement is very good for both the low temperature and 
the high temperature OPT. 

To summarize, we obtained the optimized perturbation 
theory results for the 2D and 3D LLL Ginzburg - Landau 
model in both vortex liquid and solid phases. The lead- 
ing approximant (gaussian) gives a rigorous upper bound 
on energy, while the convergent series allow one to make 
several definitive qualitative conclusions. The intersec- 
tion of the magnetization lines in only approximate not 
only in 3D, but also in 2D. The theory by Tesanovic |l^] 
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describes the physics remarkably well at very high tem- 
peratures, but deviates on the 5-10% precision level at 
ttT = — 2 in 2D and has certain imprecise qualitative fea- 
tures in the solid phase. Comparison with Monte Carlo 
simulations and some experiments shows excellent agree- 
ment. 
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FIG. 2. The 2D scaled LLL magnetization. Compari- 
son of data from Jin et al in ref. 3 with OPT calculation, 

Tesanovic et al result of ref. 10 (eq. 9) and phenomenological 
"interception" theory eq. 10 are shown for comparison. 

FIG. 3. Specific heat, 2D. Comparison of MC data 

with solid OPT (first two orders), hquid OPT (lO"* order). 
Tesanovic et. al. theory and phenomenological formula are 
also shown. 
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